Observation of non-local dielectric relaxation in glycerol 
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Since its introduction, liquid viscosity and relaxation time r have been considered to be an intrinsic 
property of the system that is essentially local in nature and therefore independent of system size. 
We perform dielectric relaxation experiments in glycerol, and find that this is the case at high 
temperature only. At low temperature, r increases with system size and becomes non-local. We 
discuss the origin of this effect in a picture based on liquid elasticity length, the length over which 
local relaxation events in a liquid interact via induced elastic waves, and find good agreement 
between experiment and theory. 



Viscosity, rj, is the basic property of matter that quan- 
tifies its resistance to flow, and was defined by Newton as 
the proportionality coefficient between shear stress and 
velocity gradient. It is commonly used to describe prop- 
erties of gases, liquids and even solids where flow is often 
referred to as "creep" [lj. Viscosity of liquids perhaps 
stands out in that list, in that theory of liquids in gen- 
eral is much less developed as compared to other phases 
of matter such as solids or gases. For example, it is stated 
2] that unlike for solids and gases, the most basic prop- 
erty of a liquid such as energy can not be calculated in 
general form. The stated reason is that interactions in 
a liquid are strong and system-specific so that the en- 
ergy strongly depends on system type. At the same time, 
atomic displacements in a liquid are large. Consequently, 
neither weakness of interactions as in a gas nor smallness 
of atomic displacements as in a solid is a small parameter 
that can be used to obtain general results 

Similarly, temperature dependence of equally basic liq- 
uid property, viscosity, has remained another major theo- 
retical hurdle in condensed matter physics. In particular, 
the behaviour of r] at low temperature is at the heart of 
the problem of glass transition [3 6]. 

A property closely related to viscosity is liquid relax- 
ation time, r. It was phenomenologically introduced 
by Maxwell in the viscoelastic picture of liquid flow as 
T) = GoqT, where is instantaneous shear modulus Q- 
Frenkel subsequently identified r with the average time 
between two consecutive atomic jumps in a liquid at one 
point in space [§]. Larger r correspond to less frequent 
atomic jumps and larger viscosities. 

Notably, in both earlier and modern [H, 0, E3- 

[Hj discussions, it was assumed that the flow property 
of a liquid, rj or r, is essentially an intrinsic property 
governed only by the liquid structure and interatomic 
interactions, and therefore is necessarily local. Applied 
to viscosity, the assertion of locality implies that as long 
as there is enough bulk liquid between two moving plates 



so that surface effects can be ignored, viscosity does not 
depend on system size. Equally, r has been considered 
to be local and independent of system size 0, 0, 0-0, 



111 . 1 1 31 ] . Consequently, several well-known relationships 
for Tj or t such as Stokes law, Einstein- Stokes relationship 
between viscosity and diffusion and so on are derived and 
discussed in the picture where rj and r are intrinsic and 
local [1, @, [IH [l3[ . The same applies to general Navier- 
Stokes equations that are used to derive and discuss a 
wide range of results and effects in liquids [HI . 

In this paper, we perform dielectric relaxation experi- 
ments in glycerol and find that r is independent on sys- 
tem size at high temperature only, but increases with 
system size at low temperature, i.e. becomes essentially 
non-local. We attribute this effect to non-local interac- 
tions between atomic jumps in a liquid via induced high- 
frequency elastic waves. We find both qualitative and 
quantitative agreement between this theory and present 
experiments. 

We first discuss the theoretical rationale for our ex- 



periment. We have recently discussed [1J, [15[ the ori- 
gin of several anomalous and intriguing relaxation laws 
in liquids that are at the centre of the problem of glass 
transition, such as Vogel-Fulcher-Tammann (VFT) law, 
stretched-exponential relaxation and dynamic crossovers. 
An important outcome of this discussion is that relax- 
ation in a liquid necessarily becomes non-local due to 
interaction between atomic jumps. The non-locality orig- 
inates as follows. 

An elementary flow process in a liquid is the jump of 
an atom from its surrounding atomic "cage" , accompa- 
nied by large-scale rearrangement of the cage atoms. We 
call this process a local relaxation event (LRE). A LRE 
lasts on the order of elementary (Debye) vibration pe- 
riod To = 0.1 ps. On this short timescale (as well as any 
other timescale that is shorter than r), the surrounding 
liquid can be viewed as an elastic medium [§]. Hence, 
large atomic motion from a LRE elastically deforms the 
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surrounding liquid, inducing elastic waves. The wave fre- 
quency, oj, is on the order of Debye frequency because the 
wavelength is on the order of interatomic separations. 
Therefore, lu > - as long as r exceeds r , i.e. in the 
whole range of t up to r = 10 2 — 10 3 s at the glass tran- 
sition temperature. As discussed by Frenkel [8j, waves 
with frequency ui > - are propagating in a liquid as in a 
solid. The waves distort cages around other LRE centres 
in the liquid, and therefore affect their relaxation. 

Hence, we identified LRE-induced elastic waves as the 
physical mechanism of mediating interactions in a liquid, 
and proposed [3] that this interaction sets the coopera- 
tivity of relaxation whose origin was widely discussed but 
not understood from the physical point of view ■ 

The key question is the range of this interaction, or 
the propagation range of high-frequency elastic waves. 
As discussed in detail |14j . this range is given by rf e i: 



del = 



(1) 



where c is the speed of sound. 

The non-trivial point is that d c \ = err increases with r. 
This is directly opposite to the commonly discussed decay 
of hydrodynamic waves, whose propagation range varies 
as 1/t. The difference is due to the solid-like regime of 
wave propagation, lot > 1, which is qualitatively differ- 
ent from the hydrodynamic regime, lot < 1 Q (see Ref. 
[14] for a detailed discussion). We called d e i liquid elas- 
ticity length because it defines the range over which two 
LREs interact with each other via induced elastic waves. 
Importantly, g? c i = cr increases on lowering the temper- 
ature because r increases. We proposed that this is the 
key to the emergence of slow relaxation, VFT law and dy- 
namic crossovers [3, EH • In this picture, the increase of 
dei and concomitant cooperativity of relaxation increases 
the activation barriers for LREs and r [LU ]. 

We note that d c \ is different from the solidity length, 
I, introduced in Ref. [l6j ]. In particular, I was derived 
not as the propagation range of high-frequency waves in 
a liquid, but as a distance that the sound wave travels 
without LREs appearing anywhere along its path during 
time -. As a result, / is proportional to (compare 
with Eq. (UJ) where c? c i oc r), is significantly smaller than 
d c i , and is limited by hundreds of nm at T g [l6| . 

As long as d e i < L, where L is system size, relaxation 
effects are not affected by system size. On the other 
hand, when d c \ > L at low temperature, one expects 
to find that relaxation depends on system size, provided 
d i < dt- Here, d t , phonon thermalization length, is the 
distance (typically up to several mm at room tempera- 
ture [17H19l |) over which an induced phonon thermalizes 
due to anharmonicity, with the result that the phonon 
spectrum of the system becomes the same as in thermal 
equilibrium. In particular, two important system size ef- 
fects are immediately predicted by this theory. First, the 



theory predicts that when d Q \ = L, r crosses over from 
the VFT law to a more Arrhenius dependence, because 
at this point further increase of d c i does not increase the 
number of interacting LREs in the system [3]. This ef- 
fect has been observed experimentally in a large number 
of glass-forming liquids [la-Hi. 

Second, this picture predicts that when d c \ > L, r in- 
creases with L and becomes essentially non-local. We 
note that in the above discussion, r can be considered 
non-local when d c \ < L as well, as long as d c i exceeds the 
distance between neighbouring LREs. However, experi- 
mental verification of the non-locality of r is most easily 
verified by varying L and demonstrating that r increases 
with L. The first demonstration of this effect in a macro- 
scopic system was done in the series of Stokes viscosity 
experiments [23} . Recently, a similar non-locality of vis- 
cosity was found in molecular dynamics simulations [3] . 
We now discuss the results of our dielectric relaxation 
experiment supporting the prediction of the size effect. 

We have measured t in glycerol, a commonly stud- 
ied glass-forming liquid, using well-known broadband 
(10-10 8 Hz in the current experiments) dielectric spec- 
troscopy technique in capacitors of two different sizes. 
Several important experimental conditions need to be 
met in order to study size effect. These include ensuring 
that temperature difference between the two capacitors 
is minimized, as well as temperature gradients in each 
system. Next, the liquid in each capacitor should be a 
closed system, and not connected to a liquid of larger 
size, as can be the case in dielectric relaxation experi- 
ments. Further, pressure variations in a closed system 
related to thermal expansion should be avoided because 
they affect r. 

To meet these requirements, we have designed and 
built the test cell shown in the inset in Figure 1. The 
cell contained two switchable capacitors CI and C2, both 
of the same radius R = 0.75 mm and different heights 
hi = 0.2 mm and = 0.02 mm, respectively. The dis- 
tance between the capacitor plates is maintained by flat 
teflon spacers. Both capacitors were cooled simultane- 
ously. Subsequently, their temperature was stabilized by 
a temperature controller with accuracy better than 0.02 
K. Due to efficient thermal coupling between the capaci- 
tors, temperature gradient through the test cell was not 
larger than 0.01 K. 

Samples of ultra-pure (>99%) glycerol were purchased 
from ICN Biomedicals. High sample purity level was ad- 
ditionally checked by refractometry. At each tempera- 
ture, the liquid was equilibrated during temperature sta- 
bilization lasting several minutes, considerably exceeding 
stress relaxation time in our temperature range. Fre- 
quency dependencies of the impedance of the capacitors 
Z(F) were measured by Quad Tech 7600 precision LCR 
Meter below I MHz and by Hewlett-Packard HP4191A 
reflectometer above I MHz. Typical dielectric loss spec- 
tra at the T =233 K, where the maximum value of 
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FIG. 1: Schematic experimental setup (inset) and a typical 
dielectric loss spectra of glycerol measured in the larger (cir- 
cles) and smaller (squares) capacitor at T=233 K. 



F m = - is in the kHz range, is shown in Figure 1. 

In Figure 2, we plot r of both capacitors as a function 
of temperature, together with log^) — log(r 2 ), where 
Ti and T2 are relaxation times in large and small capaci- 
tors, respectively. Figure 2 makes three important points 
about the predicted size effect. 

First, the size effect is absent at high temperature, 300 
K. According to Figure 2a, this temperature corresponds 
to t < 10~ 8 s and hence d cl < 0.02 mm (c f«2000 m/s for 
glycerol). Therefore, d e \ is just smaller than the smallest 
dimension of both capacitors {h-z = 0.02 mm). According 
to the theoretical prediction, r should be the same in 
both systems: t± =t%, in agreement with Figure 2b. 

Second, we observe that the size effect is present at low 



temperature, below 250 K. At 250 K, r 



10" 



(see 



Figure 2a), giving d c \ sa 2 mm. Therefore, d c \ is larger 
than the largest dimension of both capacitors below 250 
K. According to theoretical prediction, r should be larger 
in the larger system in this temperature range: t\ > 
Ti . This is in agreement with the experimental results in 
Figure 2b. 

Third and finally, we calculate the magnitude of the 
size effect from the theory. According to the theory [3] , 
the activation barrier for a LRE, U, is 



U = U + 



To fdV 

4-7T 



(2) 



The first term in Eq. (|2"]>. Uq, is the high-temperature 
intrinsic, or non-cooperative, activation barrier for a LRE 
than depends on structure and the type of interactions 
in a liquid but not on the interaction with other LREs. 




FIG. 2: (a) r in large (upper symbols) and small (lower sym- 
bols) capacitors. The horizontal long-dashed line corresponds 
to r above which the size effect is predicted to exist, (b) 
log(ri) — log(T2), where ri and T2 are relaxation times in large 
and small capacitors, respectively. 



The second term describes the effect of elastic interaction 
between LREs. Here, To is the VFT temperature, and -\ 
represents the decay of elastic stress. The integration is 
performed over the volume inside the system boundaries 
if d e \ > L an d inside d e \ if d c \ < L. We note that if 
d c i < L, U is temperature-dependent because it depends 
on d e \ and therefore on t, giving the VFT law for r [T3 |. 
The lower integration limit in Eq. @ is given by the size 
of a local rearranging region, d (do ~ 10 A). 
Integrating in cylindrical coordinates gives 



T 2hR 

U = U Q + In 

2 d {h + V/i 2 + i? 2 ) 



(3) 



where we have taken into account that do ^ R- 

We define the magnitude of the size effect, A, as 
A = — , where T\ and Ti are relaxation times in larger 
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and smaller capacitors, respectively. Using Eq. ([3]) and 



ln(r) = ln(ro) 
we write: 



where tq is Debye vibration period, 



In A = -5- In 

2T 



h x h 2 



R 2 



h 2 h! + y/h\ + R 2 



(4) 



In the low temperature range where the size effect 
is seen in Figure 2b (208-253 K), Eq. ©gives A = 
1.7 — 1.9, where we have used To = 135 K [25J and the 
above values of hi, h 2 and R. The experimental values 
of A from Figure 2b are in the range 1.2—1.3. Therefore, 
the theoretical prediction of the size effect is in reason- 
able agreement with the experimental value, providing, 
importantly, a correct order of magnitude of the effect. 
We observe that t\ and t 2 themselves vary by many or- 
ders of magnitude in the above temperature range (see 
Figure 2a), hence a reasonably good agreement of ^ be- 
tween theory and experiment is particularly encouraging. 

We note that Eq. Q predicts the decrease of A with 
temperature. Consistent with this, there is a decreasing 
trend of A in Figure 2b despite the larger error bars at 
the boundaries of the low temperature range, because 
the errors in the middle of the range were smaller in our 
experiments. 

We attribute the higher theoretical value of A to the 
decay of LRE-induced waves due to anharmonicity 
[l9l | . This decay is not considered in the theory, and leads 
to the reduction of the wave amplitude within d e \ and 
therefore weakening of the size effect. 

We finally note that T g decreases, compared to the bulk 
values, in liquids confined to nano-scale pores as well as 
in thin polymer films, implying that r becomes smaller in 
these systems [26j . This effect is therefore consistent with 
our theory. However, as reviewed in Ref. [26| in detail, 
it has remained unclear whether this is an intrinsic effect 
and to what extent the observed size effects are related 
to extrinsic factors such as the interaction between the 
liquid and the surface, particularly important in small 
systems. Our results, on the other hand, are done on es- 
sentially macroscopic samples where the above extrinsic 
effects are expected to be small. Our results therefore 
support the existence of a genuine intrinsic size effect in 
liquids, in accordance with our recent theory. 

In summary, our dielectric relaxation experiments 
show that relaxation time of glycerol increases with sys- 
tem size at low temperature, implying that r (and sim- 
ilarly, viscosity) become non-local. We propose that the 
origin of this unexpected behaviour is the long-range 
elastic interaction between local relaxation events, whose 



range increases on lowering the temperature and extends 
to system size. The results of our experiments are in a 
good qualitative and quantitative agreement with theo- 
retical predictions. 
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